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ABSTRACT 

We study the nature of the [Ti Ii] and [Ni n] emission from the so-called strontium 
filament found in the ejecta of r] Carinae . To this purpose we employ multilevel mod- 
els of the Ti ii and Ni ii systems which are used to investigate the physical condition 
of the filament and the excitation mechanisms of the observed lines. For the Ti II 
ion, for which no atomic data was previously available, we carry out ab initio calcu- 
lations of radiative transition rates and electron impact excitation rate coefficients. 
It is found that the observed spectrum is consistent with the lines being excited in a 
mostly neutral region with an electron density of the order of 10 7 cm -3 and a temper- 
ature around 6000 K. In analyzing three observations with different slit orientations 
recorded between March 2000 and November 2001 we find line ratios that change 
among various observations, in a way consistent with changes of up to an order of 
magnitude in the strength of the continuum radiation field. These changes result from 
different samplings of the extended filament, due to the different slit orientations used 
for each observation, and yield clues on the spatial extent and optical depth of the 
filament. The observed emission indicates a large Ti/Ni abundance ratio relative to 
solar abundances. It is suggested that the observed high Ti/Ni ratio in gas is caused 
by dust-gas fractionation processes and does not reflect the absolute Ti/Ni ratio in 
the ejecta of r/ Carinae. We study the condensation chemistry of Ti, Ni and Fe within 
the filament and suggest that the observed gas phase overabundance of Ti is likely the 
result of selective photo-evaporation of Ti-bearing grains. Some mechanisms for such 
a scenario are proposed. 

Key words: atomic data - atomic processes - line: formation - stars: abundances - 
stars: Eta Carinae 


1 INTRODUCTION 

r) Carinae is the most luminous star in our part of the Galaxy 
and indeed is one of the most mysterious known to us 
(Davidson & Humphreys 1997). It is surrounded by com- 
plex ejecta created in several outbursts of the star and an 
outgoing massive wind with periodic events every 5.54 years 
(Damineli 1996). The largest mass ejection occurred during 
the so-called Great Eruption observed in the 1840s - which 
the star survived. In that event, rj Carinae released about 
10 49 ' 6 ergs of kinetic energy contained in more than 10 M© 
of material (Smith et al. 2003), which has since expanded to 
form the bipolar Homunculus nebula we see today (Gaviola 
1950; Smith & Gehrz 1998; Morse et al. 2001). Despite the 
hot and very luminous central star, the circumstellar neb- 
ula is cool enough to produce a remarkably rich emission 


spectrum of lines from neutral and low ionization species 
like Fe + and Ni + (Davidson et al. 2001; Zethson 2001; Zeth- 
son et al.2001), line-of-sight absorption from several neutral 
and singly ionized atoms (Gull et al. 2005), and even emis- 
sion from H 2 and very cool dust (Smith 2002a; Smith et 
al. 2003). In addition to being dense and cool, the ejecta of 
q Carinae are nitrogen rich, representing ashes of the CNO 
cycle (Davidson et al. 1982, 1986; Smith & Morse 2004; 
Verner et al. 2005). It is the high density, low-excitation, and 
CNO abundances of the ejecta that make the Homunculus 
a unique laboratory among circumstellar nebulae, having a 
profound effect on the results discussed in this paper. 

In February 1999 a series of STIS observations of r] Cari- 
nae was made of a region about 1.5” northwest of the star 
(Zethson et al.2001). The spectrum from this region was 
significantly different from the nebular emission spectra 
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elsewhere in the Homunculus, with two of the lines being 
identified as forbidden emission of singly-ionized strontium 
([Sr II]). The measured velocity of the filament (~ —100 
km/s) and its distance from the central star are consistent 
with the gas having been ejected from the same star during 
the Great Eruption. 

Further spectroscopic observations and a complete list 
of lines identified in the filament were reported by Hartman 
et al. (2004). They show that lines of Fe I, V II, and Ti II 
are strong in the filament, whereas only a few Fe II lines 
are detected. Other lines in the spectrum at about the same 
velocity come from C I, Mg I, Al II, Ca I, Ca II, Sc II, Gr II, 
Mn II, Co II, and Ni II. Curiously no H, He, S, O, or N lines 
are present in this region. 

The cosmic abundance of strontium is lower than that 
of iron and nickel by factors of 30,000 and 3,000 respec- 
tively. Thus, strontium emission is normally undetectable in 
nebular spectra. The fact that [Sr II] emission was detected 
with strengths comparable to those of [Fe il] and [Ni II] lines 
indicates that either strontium is locally overabundant in 
r i Carinae by orders of magnitude, which is inconsistent with 
the evidence of CNO-cycle-processes ejecta, or very unusual 
conditions exist in this region of the nebula. 

In a recent study of the Sr II and [Sr n] spectra of the 
filament (Bautista et al.2002) we found that [Sr n] emission 
may be strongly enhanced by continuum fluorescence, which 
partially explains the strength of these lines. Although it was 
not possible at that time to establish the precise abundance 
of Sr, this result did imply that Sr might not be as over- 
abundant as previously thought. This point opened up new 
questions. If [Sr II] lines can be efficiently excited by fluores- 
cence, why are the lines seen only in this structure, but not 
elsewhere in the Homunculus, and not in any other emission 
line nebulae? 

The observed spatial and velocity distribution of the Sr- 
filament suggests that it is located in the equatorial region 
between the hourglass-shaped bipolar Homunculus (David- 
son et al 2001: Smith 2002a; Ishibashi et al, 2003). The Sr- 
filament is located at a projected distance of 10 lightdays 
(at 2300 pc distance), where it should receive intense UV 
radiation from 77 Carinae . The spectra of the Weigelt Blobs 
B and D, located a few light days from 77 Carinae , require 
the spectral flux distribution of an O star with Te=35,000K 
during the -5-year broad maximum (Verner et al.2005). The 
Sr-filament should receive a diluted amount of this same UV 
radiation, but its emission line spectrum indicates a much 
lower ionization condition. The Sr-filament appears to be 
shielded not only by neutral hydrogen (13.6 eV), but also 
singly-ionized iron (7.9 eV) which requires a substantial col- 
umn density of Fe 11 between 77 Carinae and the Sr-filament. 
One possibility is that the Sr-filament is located at the outer 
side of an apparent optically thick, dusty, torus seen in the 
mid-IR (Smith et al. 2003). In order to demonstrate this 
scenario, we must determine the nature of a hot star contin- 
uum filtered out by an Fe II blocking column and predict the 
physical conditions within the filament that would sustain 
such a low ionization region. Then we can determine the 
actual abundances of all species observed in this structure. 

In this paper we study the emission spectra of two 
prominent ions in the filament: Ti II and Ni II. Ti II is stud- 
ied for the very first time as the relevant atomic data has 
never before been obtained. Ti 11 emission lines are rarely 


Table 1. Ti II configurations and scaling parameters 


Spectroscopic configurations 

3s 2 3p 6 3d 3 , 3.s 2 3 p 6 3d 2 4s, 3s 2 3p 6 3d 2 4 p, 3s 2 3p 8 3d 4s 2 
Correlation configurations 

3s 2 3p 6 3d 2 id, 3s 2 3 p 6 3d 4s 4p L 3.s 2 3p 6 3d 4s 4d,_ 

3s 2 3p 6 3d 4s 5s, 3s 2 3p 6 3d 4s 5d, 3s 2 3 p 6 3d 4s 5p, 

3s 2 3p 6 3d 4s if, 3 s 2 3p 6 3d 4p 2 , 3s 2 3p 5 3d 4 , 3s 2 3p 5 3d 3 4s, 
3s 2 3p 5 3d 3 4p, 3s 2 3p 5 3d 2 4s 2 , 3s 2 3p 4 3d 5 , 3s 2 3p 4 3d 3 4s 2 
3s 2 3p 4 3d 4 4s, 3s 2 3p 4 3d 4 ip, 3s 3 p 6 3d 4 , 3s 3p 6 3d 3 4s, 

3s 3p® 3d 3 4p 3s 3 p® 3d 4 4s, 3s 3 p 5 3d 4 4p, 3p 6 3d 4 4s, 

3p 6 3d 4 4p 


Ail 

1.43570(ls), 1.1 2290 (2s), 1.06650(2p), 1.09690(3s), 1.08430(3p), 
1.12590(3d), 1.27470(4s), 1.22660(4p), -0.79940(4_d) 

-0.64020(5s) , 2.99520(5d), 1.90030(5p), 1.88930(4/) 


seen in emission nebulae as its ionization potential, 13.58 
eV, leads to Ti HI in most H II regions. Ni 11 , often seen in 
nebular emission spectroscopy, is much better understood 
and herein serves as a reference ion. 

The remainder of the paper is organized as follows: In 
Section 2 we present the atomic models built for the Ni 11 
and Ti 11 systems, with particular detail in the calculations 
of atomic parameters that had to be performed for Ti II. In 
Section 3 we show the results of the nebular diagnostics. In 
Section 4 we derive the relative abundances of Ti 11 and Ni 11 
in gas phase, while in Section 5 we discuss possible scenarios 
of Ni and Ti depletion in the nebula. Finally, our discussion 
and conclusions are given in Section 6 . 


2 ATOMIC DATA AND ATOMIC MODELS 
FOR Ni n AND Ti 11 

For the present analysis of [Ni 11 ] spectra we employ a 76- 
level model for the ion that includes continuum fluorescence 
excitation as in Bautista et al.(1996) and improved colli- 
sions! data reported by Bautista (2004). This new collisional 
data, that results from very extensive calculations, is ex- 
pected to be more accurate than previous data. 

2.1 Radiative data for Ti 11 

For Ti II, neither complete transitions rates for forbidden 
transitions nor collision strengths had been reported, thus 
large calculations of atomic data for this ion were carried 
out within the present work. 

For the calculation of radiative transition probabilities 
and for the wavefunctions of the Ti II system we use the 
atomic structure code AUTOSTRUCTURE (Badnell 1986) 
to reproduce the structure of the ion. This code is based on 
the program SUPERSTRUCTURE originally developed by 
Eissner et al. (1974), In this approach the wavefunctions are 
written as a configuration interaction expansion of the type: 

(i) 

3 

where the coefficients Cji are chosen so as to diagonalize 
(ipi | H | where H is the Hamiltonian and the ba- 
sic functions 4>j are constructed from one-electron orbitals 
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Table 2. Energy terms in Rydberg for the Ti II ion relative to the 
ground state. The table compares experimental energies with ab 
initio calculated energies (Theo.) and energies with term energy 
corrections (Theo .-TEC) 



Config. 

Terra 

Theo. 

Theo.-TEC 

Exp. 

1 

3d 2 { 3 F)4s 

a 4 F 

0.0 

0.00 

0.0 

2 

3d 3 

b 4 P 

0.08806 

0.00834 

0.007839 

3 

3d 2 ( 3 F)4s 

a 2 F 

0.06439 

0.04231 

0.041528 

4 

3d 2 { 1 D)4s 

a 2 D 

0.08922 

0.05534 

0.077508 

5 

3d 3 

a 2 G 

0.09982 

0.08140 

0.080552 

6 

3 d 3 

a 4 P 

0.09202 

0.08543 

0.084077 

7 

3d 3 

a 2 P 

0.09217 

0.07588 

0.088476 

8 

3d 2 ( 3 P)4s 

b 4 P 

0.10519 

0.09021 

0.088784 

9 

3d 3 

b 2 D 

0.13436 

0.11686 

0.113737 

10 

3d 3 

a 2 H 

0.14976 

0.11788 

0.113935 

11 

3d 2 ( 1 G)4s 

b 2 G 

0.18304 

0.14025 

0.137017 

12 

3d 2 l 3 P)4s 

b 2 P 

0.19358 

0.16370 

0.149115 

13 

3d 3 

b 2 F 

0.22551 

0.19290 

0.188561 

14 

3d4s 2 

c 2 Z) 

0.37419 

0.20894 

0.226677 

15 

3d 2 ( 3 F)4p 

2 4 G° 

0.24771 

0.26968 

0.270750 

16 

3d 2 ( 3 F)4p 

2 4 F° 

0.25782 

0.28059 

0.281422 

17 

3d 2 { 3 F)4p 

z 2 po 

0.26630 

0.27996 

0.285498 

19 

3d. 2 ( l S)4s 

a 2 S 

0.35738 

0.28071 

0.287619 

18 

3d 2 ( 3 F)4p 

z 2 D° 

0.27074 

0.29066 

0.288805 

20 

3d 3 

d 2 D1 

0.28882 

0.34006 

0.292376 

21 

3d 2 ( 3 F)4p 

z 4 D° 

0.27443 

0.29206 

0.295841 

22 

3d 2 ( 3 F)4p 

z 2 G° 

0.30665 

0.31108 

0.313768 

23 

3d 2 { 3 P)4p 

z 2 S° 

0.32376 

0.33886 

0.339041 

24 

3d 2 { l D)4p 

z 2 P° 

0.34940 

0.35509 

0.356748 

25 

3d 2 ( 1 D)4p 

y 2 D° 

0.35800 

0.35671 

0.358146 

26 

3d 2 ( 1 D)4p 

y 2 F° 

0.35890 

0.35852 

0.362557 

27 

3d 2 ( 3 P)4p 

z 4 S° 

0.35195 

0.36370 

0.362702 

28 

3d 2 { 3 P)4p 

y 4 D° 

0.36084 

0.36767 

0.368031 

29 

3d 2 l 3 P)4p 

z 4 P° 

0.37814 

0.38047 

0.381834 

30 

3d 2 ( 1 G)4p 

y 2 G° 

0.40620 

0.39781 

0.396745 

31 

3d 2 ( 3 P)4p 

X 2 D° 

0.41916 

0.41092 

0.407173 

32 

3d 2 ( 3 P)4p 

y 2 P° 

0.42858 

0.40434 

0.412787 

33 

3d 2 { l G)4p 

z 2 H° 

0.43455 

0.41952 

0.415324 

34 

3d 2 ( 1 G)4p 

x 2 F° 

0.45445 

0.43599 

0.431113 

35 

3d 2 { 1 S)4p 

2 po 

0.61270 

0.57783 

0.569851 


generated using the Thomas-Fermi-Dirac model potential 
(Eissner & Nussbaumer 1969). The A n i scaling parameters 
in this model are optimized by minimizing a weighted sum 
of energies. Representing a system as complex as Ti II re- 
quires many configurations, which must be selected rather 
carefully as to account for all strongly coupled configura- 
tions without exceeding available computational resources. 
Experience shows that configurations with one and two elec- 
tron promotions out of the 3s, 3p, and 3d orbitals are often 
important. Further, improvements can be obtained in some 
cases through the use high laying non-physical orbitals. For 
the present Ti II model it was necessary to include several 
non-physical orbitals in the expansion, i.e. 4d, 5 id, 5s, 5p, 
4/, The configuration representation and scaling parame- 
ters used in this work are listed in Table 1. 

Fine tuning (semi-empirical corrections) — which is use- 
ful for treating states that decay through weak relativis- 
tic couplings (e.g. intercombination transitions) — takes the 
form of term energy corrections (TEC). By considering the 
relativistic wavefuntion, ipl, in a perturbation expansion of 
the non- relativistic functions 


= ipf r + ^2 4>j 


j^i 


(^ nr |gib + gabion 

E nT - E nl 


(2) 


where H ib and H 2 b are the one- and two-body Breit-Pauli 
relativistic operators. Thus, a modified H nr is constructed 
with improved estimates of the differences Ef r — E™ so as 
to adjust the centers of gravity of the spectroscopic terms 
to the experimental values. Table 2 shows a comparison be- 
tween the observed energies averaged over the fine struc- 
ture taken from Huldt et al. (1982), Moore (1949), and Rus- 
sell (1927) and compiled in the NIST energy levels database 
(NIST 2000), the computed ab initio term energies, and the 
computed energies after TEC corrections. The overall agree- 
ment between ab initio energies and experimental values is 
typically within 10%, with better agreement for the odd par- 
ity terms but with some problems for low laying even parity 
multiplets. To the limit of our computations, we were unable 
to find a set of correlation configurations able to reproduce 
ab initio the correct order of the terms 3d 3 a 2 G, a 4 P, 
and a 2 P (terms 5 through 7). We also found it difficult 
to reproduce the position and an accurate energy for the 
3d4s 2 c 2 D term. These large discrepancies between calcu- 
lated and experimental energies were readily resolved by the 
TEC technique, although some corrections were uncomfort- 
ably large. A noticeable exception is the a 2 D term which 
lowers its energy far from the experimental value after TEC 
corrections. This term could not be improved further as it 
has a profound effect on the perturbative corrections on the 
b 4 F and a 2 F terms. 

This ion’s representation was used to compute radia- 
tive transition rates for both dipole-allowed and forbidden 
transitions. Table 3 shows a comparison between the present 
log(gf)-vahies for dipole allowed transitions and those de- 
termined experimentally by Pickering et al. (2001, 2002) 
and Bizzarri et al. (1993), computed by Kurucz (2000), and 
those compiled in the NIST database (taken from Meggers 
et al. 1975; Roberts et al. 1975; Danzmann and Kock 1980; 
Blackwell et al. 1982; Roberts et al. 1973). Most gf-v , alues 
in the NIST database have an accuracy rating ”D”, which 
means that the uncertainties are greater than 50%. Ta- 
ble 3 is limited to lines with <? /-values greater than 0.01 for 
which reasonably accurate data is expected given the present 
representation of the ion. We find good agreement, within 
0.2 dex, between our results and those of other authors for 
log{gf) > — 1 and a mean dispersion of about 0.5 dex for 
smaller values. Given the uncertainties in the present re- 
sults, for our Ti II multilevel model we combine the NIST 
/-values whenever available with the present radiative data 
for all other transitions. 

Our representation of the ion was also used to compute 
electric quadrupole and magnetic dipole transition rates. 
The results of such calculations are rather sensitive to con- 
figuration interactions and level mixing. Some indication 
about the accuracy of the present rates is obtained from 
recent measurements by the FERRUM project (Hartman et 
al. 2003, 2005) of lifetimes of metastable levels, which de- 
cay through parity forbidden transitions. Lifetimes for five 
levels have been reported so far, which are compared in Ta- 
ble 4 with our theoretical results. Our calculated lifetimes 
agree within ~50% with the experimental determinations, 
which is the estimated overall uncertainties for the present 
radiative rates for forbidden transitions. 
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Table 3. Comparison between log(g f)-vaXues for dipole allowed 
transitions from the present work, and those of Pickering et al. 
(2001, 2002), Bizzarri et al. (1993), Kurucz (2000), and the values 
recommended by NIST 


Wave, (nm) 

Present 

Pickering 

Bizzarri 

Kurucz 

NIST 

Upper level: 29544 cm 1 

3d 2 ( 3 F)4p 

z 4< J f s/2 



401.23835 

-1.14 

-1.84 

-1.75 

-1.73 

-1.61 

349.10494 

-1.01 

-1.15 

-1.06 

-1.13 


340.98043 

-1.80 

-1.98 

-1.89 

-1.89 

-1.90 

339.45720 

-0.56 

-0.55 

-0.54 

-0.52 

-0.59 

338.37587 

-0.052 

0.16 

0.15 

0.20 

0.050 

Upper level: 31490 cm 1 

3d 2 ( 3 F)4p 

2 2jC V/2 



465.72004 

-2.68 

-1.59 

-1.63 

-1.59 

-2.15 

439.50311 

-0.92 

-0.54 . 

-0.51 

-0.45 

-0.66 

375.92914 

0.16 

0.28 

0.27 

0.31 

0.20 

321.47669 

-3.30 

-1.40 

-1.34 

-1.35 

■ -1.39 

319.75185 

-0.30 

-1.97 

-1.90 

-2.55 

-1.87 

1587.3644 

-1.38 



-1.92 



Table 4. Lifetimes (in sec) for metastable levels of Ti H as deter- 
mined experimentally by Hartman et al. (2003, 2005) and from 
the present calculations 


Level 

Hartman 

present 

3d 2 ( 3 P)4s b 4 F 3/ o 

18(4) 

31.1 

3d 2 ( 3 P)4sb i P 5/2 

28(10) 

16.4 

U 2 ( 3 P)4s b 2 P 1/2 

14(3) 

12.7 

3d4s 2 c 2 Av 2 

0.29(0.01) 

0.192 

3d4s 2 c 2 D 6/2 

0.33(0.02) 

0.196 


2.2 Collisional data for Ti 11 

We have computed electron impact collision strengths for 
transitions among the lowest 82 energy levels of Ti II that 
belong to the configurations 3d 3 , 3d 2 4s, 3d 2 4p, and 3d4s 2 . 
We employ the Breit-Pauli R-matrix method and the BPRM 
set of codes developed in the framework of the IRON Project 
(Hummer et al. 1993; Berrington et al. 1995). Partial wave 
contributions are included from 60 SLrr total symmetries 
with angular momentum L = 0 — 9, total multiplicities 
(25 + 1) = 1 — 5, and parities even and odd. Then, the colli- 
sion strengths were ’’topped up” with estimates of contribu- 
tions from higher partial waves based on the Coulomb-Bethe 
approximation (Burgess 1974). Collision strengths for the 
fine-structure levels were obtained by algebraic recoupling 
of the LS reactance matrices (Hummer et al.1993). The col- 
lision strengths were calculated at 11000 energy points from 
0 to 6 Ry, with 90% of the points in the region with res- 
onances. This number of energy points was found to pro- 
vide sufficiently good resolution for accurate calculations 
of the contributions of autoionizing resonances to the ef- 
fective collision strengths at temperatures between 5000 K 
and 20000 K. Figure 1 shows the collision strengths for some 
of the dominant transitions. 

We compute Maxwellian-averaged effective collision 
strengths, defined as: 


T(T)= / Qif(e/)exp(—€f/kT) d(e//kT), (3) 

Jo 

where Q is the collision strength for the transition i to 
/, and «/ is the energy of the outgoing electron. The T’s 
were computed for various temperatures between 5000 and 
20 000K to be used in the calculation of excitation rates. 

2.3 Multilevel model of Ti II 

We built an 82-level excitation equilibrium model for Ti II 
which considers both electron impact excitation and contin- 
uum fluorescence excitation. A model of this size is large 
enough to account for collisionally excited emission at the 
temperatures of interest (~ 10 4 K), and the dominant pho- 
toexcitation channels from the ground level. If more levels 
were included in our models they would tend to increase the 
efficiency' of continuum fluorescence excitation, although it 
has been shown that the contribution of individual levels de- 
creases rapidly with the level energy (Bautista et al. 1996). 
We assume that the radiation field density, U v , at photon 
energies below the Lyman ionization limit (13.6 eV) can be 
approximated by a blackbody with temperature Tr times a 
geometrical dilution factor w, i.e. 

c 3 U v _ w 

S-n-hu 3 sxp{hu/kTR) — 1 

For the present calculations we adopt a blackbody temper- 
ature of 35 OOOK (see Verner et al.2005) while the dilution 
factor is diagnosed from the observed spectra. A different 
choice in the blackbody temperature would lead to a change 
in the diagnosed dilution factor, but would have negligible 
effect on the modeled spectrum. 

Fig. 2 presents a diagram of the dominant photo- 
excitation channels for two [Ti IlJ lines arising from quar- 
tets and doublets terms. The main photoexcitation chan- 
nels are indicated by dashed lines and the radiative deex- 
citation transitions are marked by solid lines. Such dipole 
allowed radiative decays lead to emission lines, observable 
in the spectra, just as they populate the levels that give 
rise to forbidden lines. It is worth noticing that in the ab- 
sence of collisional excitation, the strength of forbidden lines 
would be tightly bound to the intensity of the dipole transi- 
tions from the pumping mechanism. On the other hand, the 
contribution of collisional excitations to the forbidden lines 
introduces a complex temperature and density dependence 
in the relative strengths of dipole and forbidden transitions, 
which can be used for spectral diagnostics as illustrated in 
the next section. 

Fig. 3 presents the line emissivity per ion of 
the Ti II lines [A6648] (b 2 G 7 / 2 ~ a 4 F 7/ / 2 ) and 

[A4918] (c 2 D 3 / 2 — a 2 F 5 / 2 ) and the blend [A100696] 
(6 4 P 3 /2 - a 4 F 3/2 + b 4 P 5/2 - a i F 5/ 2 ) vs. N e for an 
electron temperature of 6000 K, a blackbody temperature 
of 35 000 K and four values of the dilution factor, w = 
0, 10 — 10 , 5 x 10 -10 , and 10~ 9 . For electron densities less 
than N e ~ 10 7 cm~ 3 continuum fluorescence dominates the 
strength of most tines. On the other hand, for electron densi- 
ties around N e ~ 10 7 , as previously diagnosed for the stron- 
tium filament (Bautista et al. 2002), collisional excitation 
dominates the [A6648] and [A100696] features, while fluores- 
cence still makes a very significant contribution to the exci- 
tation of the c 2 Z? 3/2 level, responsible for the [A4918] line. 
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0 0.2 0.4 0.6 0 0.2 0.4 0.6 

Energy (Ryd) 


Figure 1 . Collision strengths vs. incident electron energy for excitations from the ground level a 4 F3/2 to a 4 F 3 / 2 (a), a 4 F 7 / 2 (b), 
« 4 -Fb/2 (c), b 4 F 3/2 (d), a 2 F 7/2 (d), a 2 D 3 /2 (e). 

This is because collisional transition rates from the ground tics were reported in Bautista et al. (2002) from the [Sr 11] 
levels to the c 2 D multiplet are rather small. lines. There we estimated an electron density of the order 

of 10 7 cm -3 and an upper limit to the dilution factor for 
photoexcitation of ~ 10~ 9 . 

3 SPECTRAL DIAGNOSTICS 

We use the [Ni II] and [Ti ll] emission line strengths to diag r 3 x observed line intensities 
nose the physical conditions of the filament on each of the 

three HST/STIS observations: March 2000, April 2001, and For the present work, we use the deep spectra extending 
November 2001 (hereafter marOO, aprOl, novOl respectively; from 1640A to 10100A analyzed by Hartman et al. (2004). 

Bautista et al.2002; Hartman et al.2004). Previous diagnos- The original data was recorded during three visits with the 
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Figure 2. Dominant photoexcitation channels for two [Ti tl] 
lines. Some of the dipole allowed lines also emitted in the photo- 
pumping process are also indicated 


Hubble Space Telescope/Space Telescope Imaging Spectro- 
graph (HST/STIS) under several different observational pro- 
grams. For the marOO observation the long 52”x0.2” aper- 
ture was aligned close to the axis of the Homunculus (Po- 
sition Angle -41deg), whereas the aperture was aligned at 
22 deg in aprOl and -130 deg in novOl. The different an- 
gles were defined primarily by HST sun-angle restrictions, 
but were used to study the extent of the Strontium filament 
in different directions. The intersection of these apertures 
(see Fig. 1 in Hartman et al.2004) is on the apparent peak 
of the Strontium filament as first noticed by the discovery 
observations (Zethson et al.2001). 

The line fluxes were integrated across the line profiles of 
the components centered at ~-100 km/s at specific positions 
along the slit. The intersection position of the different visits 
is the primary position of interest. However as noted below, 
several additional positions were sampled along the aperture 
oriented radially to rj Carinae (P.A.=-41°) and tangentially 
(-130°) to understand the spatial variation of specific line 
ratios and therein, to study the excitation and ionization 
structure of the filament. The measured line intensities vary 
along the aperture and between visits. The structure is com- 
plex as other nebular structures are in line of sight: namely 
the Homunculus and thfe Little Homunculus. At some posi- 
tions, Hartman et al.(2004) noted a second velocity structure 
in many lines. These two velocity systems are likely related 
to the two disks noted by Davidson et al.(2001) located be- 
tween the bipolar lobes. The two velocity components cor- 



Figure 3. Line emissivity of the [Ti uj A6648, A10069 and A4918 
lines ion vs. electron density for dilution factors, w = 0 (solid 
line), 10 -10 (dotted line), 5xl0~ l ° (short-dashed line), and 10 -9 
(long-dashed line). The curves are calculated for T e = 6000 K, 
N e = 10' cm -3 , and a blackbody temperature of 35 000 K. 


respond to material in the same equatorial disk but with 
two different ages, as noted by Davidson et al. (2001). The 
more prominent feature is associated with the parts of the 
skirt ejected along with the Homunculus in the 1840’s. The 
fainter and more blueshifted feature is also in the equator, 
but was ejected in the 1890 outburst along with the Lit- 
tle Homunculus (Ishibashi et al. 2003; Smith 2005) and the 
Weigelt knots (Smith et al. 2004; Weigelt et al. 1995). In 
this study we focus on the structure ejected in the 1840s. 

As with any spectroscopic analysis the observed line in- 
tensities must be corrected for extinction. None of the lines 
that form the classical line ratios (e.g. hydrogen recombina- 
tion lines) are detected. We use two relatively well known 
[Ni il] transitions that arise from a common upper level, i.e. 
A7413(a 2 £> 3/2 - a 2 F 5/2 ) and A6668(a 2 D 5 / 2 — a 2 -£5/2)1 
and the standard interstellar extinction curve of Cardelli 
et al.(1989). While these lines are relatively close in wave- 
length, thus far they are the best candidates for estimat- 
ing extinction as all other lines of interest are much weaker 
or their predicted line ratios are less certain. The observed 
line ratios are 2.20 ± 0.13 in marOO, 1.94 ± 0.13 in aprOl, 
and 2.07 ± 0.13 in novOl. Thus, there are no significant 
variations of the line ratio between the three observations. 
From these observations we derive an extinction magnitude, 
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Av = 2.1 ± 0.9, as required to correct the ratios to the 
theoretical value of 1.6, with an estimated uncertainty of 
~10%. The relatively large Uncertainty in Av comes from 
the fact that the two lines used in the diagnostic are close 
in wavelength. This value for the extinction magnitude is 
consistent with the value measured, Av = 2.06, from pho- 
tometry of stars of the Tr 16 cluster in the Carinae nebula 
(Tapia et al.2003). The actual extinction magnitude to the 
filament could be somewhat higher than to the Carinae stars 
due to the dust within the nebula. Smith (2002b) studied lo- 
cal extinction within the larger Carina Nebula, and found 
a total-to-selective extinction of R—Av/ E(B — V) = 4.8, 
which is different than the normal ISM value of 3.1. Fur- 
thermore, the extinction around r] Carinae itself is likely to 
be very gray because of the large grain size (Smith et al. 
2003). However, thermal-IR maps of the dust column den- 
sity in the Homunculus reveal a local minimum near the 
position of the Sr filament (Smith et al. 2003), so the un- 
usual reddening there may not be a severe problem. We thus 
adopt for the remaining of this paper a extinction law with 
Av = 2.1 and the interstellar R = 3.1. Clearly, the remain- 
ing uncertainty on this value can be important in trying to 
compare lines well separated in wavelength, thus all diag- 
nostics should be carried out with lines close to each other. 

Tables 5 and 6 present the complete lists of [Ni n] and 
[Ti n] lines used for the present study and their intensities 
corrected for extinction. 



N < =10 a cm* 





3.2 [Ni it] diagnostics 

The combined contributions of collisional excitation and 
fluorescence make spectral diagnostics much more difficult 
than with traditional collisionally excited lines. Most line ra- 
tios of interest become simultaneously sensitive to electron 
density, temperature, and the continuum radiation field (i.e. 
dilution factor), thus many lines have to be studied as a 
whole to derive reliable conclusions about the physical con- 
ditions of the emitting plasma. 

In the case of [Ni II] lines, only a few diagnos- 
tic line ratios can be used because out of the seven 
dipole forbidden lines measured in the spectrum sev- 
eral share the same upper level. Fig. 4 shows the 
line ratio /([AA6668 + 7413])//([A7379]), /([ AA4327 + 
7257] ) //([A7379] ) , and /([AA4327 + 7257])//([AA6668 + 
7413]) against the logarithm of the dilution factor. The ra- 
tios were computed for three values of the electron density, 
N e = 10®, 10 7 , and 10 s cm ~ 3 displayed in separated panels 
from left to right in the figure. In each panel we present the 
line ratios as computed for temperatures of 5000, 7000, and 
9000 K. The observed lines ratios are indicated in the plots 
by horizontal lines for each observation. Two conclusions can 
be derived from these plots in regards to the dilution factor 
in the filament. First, for the novOl observation the value 
of logio w is in between -9.2 and -8.5. There is evidence for 
different dilution factors for each of the three observations, 
with a difference as large as 1 dex between the marOO and 
the aprOl observations, with the higher value correspond- 
ing to the marOO observation. The dilution factor from the 
novOl observation lies between the values for the other two 
observations. The spread in line ratios, and consequently on 
derived physical conditions, in the three observations are 
the result of different orientations of the spectrograph slit 


Figure 4. [Ni ll] emissivity line ratios 
/([A A6668 +7413] ) // ( [A7379]) , 7([AA4327+7257])//([A7379]), and 
/([AA4327+7257])//([AA6668+7413]) against the logarithm of the 
dilution factor. The line ratios computed at N e = 10®, N e = 10' , 
N e = 10 8 cm -3 are presented in separate panels located from 
left to right in this figure. For each value of N e the line ratios are 
computed for temperatures of 5000, 7000, and 9000 K. The ob- 
served line ratios are indicted by horizontal lines as dashed-line: 
marOO; dotted- line : aprOl; solid line: novOl 


with respect to the star (see Hartman et al. 2004). For the 
marOO observations the slit was centered on rj Carinae with 
the fiducial, F2, blocking the starlight. This provided a spa- 
tial sampling of the Sr-emitting region extending radially 
from the central source. That this spectrum yields the high- 
est value for the dilution factor is due to finding the ra- 
dial distance for which both apply. The novOl spectrum was 
taken at position angle -130 degrees and provided a slice of 
the Sr-filament that is tangential to the distance from the 
central source. As the position of the slit is more distant 
than the optimal position for highest dilution factor, we get 
a lower value for the dilution factor along this slit. The di- 
lution factor changes by an order of magnitude between the 
optimal radial position and the tangential sample because 
these observations sample slightly different spatial regions. 
Photons with energies below 7.9 eV encounter increasing 
optical depths with distance from the central source. For 
these energies the main source of opacity is free-free opacity 
(Bremsstrahlung) . Assuming an electron temperature of the 
order of 10 4 K, density ~ 10 7 cm -3 , and a difference in the 
path length within the filament between the two observa- 
tions of ~ 0”.3 « 10 16 cm, one gets an optical depth, r 
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Table 5. [Ni II] lines from the Sr-filament 


lab 

Transition 

Observed flux (10 13 
ergs sec " 1 cm ” 2 arcsec -1 ) 

Extinction corrected flux 

(Ay = 2 . 1 ) 




MarOO 

AprOl 

NovOl 

MarOO 

AprOl 

NovOl 

7379.970 

8303.270 

a 2 F 7 / 2 
a 2 1*7/2 

— a ^£^ 1/2 
— a 2 £> 3/2 

286.00 

211.00 

239.00 

11.80 

1088.0 

802.7 

909.2 

34.6 

6668.670 . 
7413.690 

a 2 1*5/2 
a 2 f*5/2 

— a 2 £^i /2 

— a 2 £> 3/ 2 

64.60 

142.0 

41.60 

80.60 

51.30 

106.0 

.304.3 

534.8 

195.9 

303.5 

241.6 

399.2 

4327.453 

7257.860 

a 4 B 5 /2 

g 4 Bo/2 

— a ^&i /2 
-a 4 F 7/ 2 

23.00 

18.10 

12.80 

27.00 

10.70 

90.8 

246.8 

50.5 

368.1 

42.2 

6815.450 

“ ' l Py2 

- a **F 5/2 

15.00 

5.71 

7.32 

67.6 

25.7 

33.0 

6793.350 

a 4 F l/2 

- a 4 F 5/2 

10.40 

3.50 

19.70 

47.2 

15.9 

89.4 


« 2, which is consistent with the observed variation in the 
dilution factor. 

3.3 [Ti n] diagnostics 

The large number of Ti 11 and [Ti n] lines present in the 
spectra allows us to constrain the physical conditions of the 
filament 

As discussed earlier, the combined contributions of col- 
lisional excitation and fluorescence lead to complex temper- 
ature and density dependence in the relative strengths of 
dipole to forbidden emission lines. Nonetheless, such ratios 
are much more sensitive to the strength of the dilution factor 
than to other parameters, so they can be used to diagnose 
the value of the dilution factor. 

In Fig. 5 we plot Ti II line ratios between dipole allowed 
and forbidden lines against the dilution factor, together with 
the observed ratios. These line ratios are very sensitive to 
the value of the dilution factor, but depend relatively little 
on T e and N e . From this figure, the line ratios for the novOl 
observation are consistent with logio w between -10.0 and 
-8.6, while the ratios from the marOO spectrum indicate a 
higher dilution factor, by as much as 0.5 dex. These results 
are in agreement with those obtained with the [Ni il] lines. 

Considering all possible line ratios among forbidden 
lines we found one that is primarily sensitive to the dilution 
factor, and depends little on the temperature and density. 
This ratio is shown in Fig. 6, as computed for T e = 5000, 
6000, and 7000 K and N e = 10 a , 10 7 , and 10 s cm' 3 . Once 
again the observed line ratios lay in the range of logio w 
between -9.1 and -8.6. 

From the agreement between various Ni II and Ti II 
diagnostics we find the value of the dilution factor to be 
log io w = -9.0 ± 0.5. With the dilution factor constrained, 
one can diagnose other parameters of the emitting region. 
Fig. 7 presents line ratios among Ti II forbidden lines vs. 
the electron temperature. The line ratios are shown in var- 
ious panels for N e = 10 6 , 10 7 , and 10 s cm -3 . In each case 
we show the line ratios for log io w = —9.0, -8.5, and - 
9.5. The observed line ratios from the novOl spectrum are 
also indicated. Inspection of these plot reveals that the 


best agreement for all the lines observed is obtained for 
T, = 6000 ± 1000 K. 

This relatively low temperature is consistent with the 
observed ionization structure of Fe and Ni in the filament. 
The observed iron emission arises mostly from neutral iron, 
indicating that this dominates over Fe II. Because of the 
iron shielding of the filament it is expected that there would 
be few photons capable of ionizing Fe i. But, even in the 
absence of photoionization collisional ionization drives most 
of the iron into Fe n whenever electron temperatures ex- 
ceed ~ 7000 K (see Arnaud and Raymond 1992). Below this 
temperature neutral iron may be the dominant ionization 
stage. On the other hand, under collisional ionization con- 
ditions the fraction of ionized to neutral nickel is more than 
twice the fraction for iron at the same temperature. Thus, 
Ni II is expected to remain the dominant ionization stage for 
temperatures as low as ~ 5000 K. 

Finally, we can determine the electron density. Fig. 8 
depicts four line ratios as function of the electron density for 
a fixed logio w=-9.0 and temperatures of 6000±1000 K. The 
electron density is readily diagnosed at around 10 7 cm -3 . 

Given the physical conditions diagnosed above we com- 
pare the overall agreement between the computed spectrum 
and observations. This comparison for the novOl spectrum 
is shown in Table 7. The agreement is reasonably good con- 
sidering the remaining uncertainties in the atomic data and 
the shape of the pumping radiation continuum. 

3.4 Spatial variations of the dilution factor 

The Ti II f([AA6652 + 6727j)//([A4918]) ratio is large across 
most of the spectral slit, which allows us to further inves- 
tigate the spatial variations of the dilution factor near the 
Sr-filament. To do this we integrated the intensity over 0.13" 
wide regions of the spectrum from marOO and novOl. The 
measured line ratios, corrected for extinction, and the de- 
rived dilution factors from them are shown in Fig. 9. The 
marOO observations, when the F2 fiducial (0.85” long) was 
centered on 77 Carinae with the large aperture extending ra- 
dially from the star at the middle of the filament (1.3” from 
77 Carinae ), clearly show a linear decay of the dilution factor 
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Table 6. [Ti It] lines from the Sr-filament 


\lab 

Transition 

Observed flux (10 13 
ergs sec -1 cm -2 arcsec -1 ) 

Extinction corrected flux 
(A v = 2.1) 



MarOO 

AprOl 

NovOl 

MarOO 

AprOl 

NovOl. 

4918.247 

c 2 -°3/2 ~ a 2 F5/2 

13.60 

8.62 

11.50 

124.8 

61.5 

105.6 

4984.205 

c 2 ^3/2 ~ a 2 FV/2 

2.64 

1.32 

2.00 

23.3 

11.7 

17.7 

6153.623 

c 2 ^3/2 ~ a 2 Dz/2 



3.76 



24.0 

4927.288 

c 2 ^5/2 ~ a 2 fy/ 2 

14.70 

13.2b 

13.40 

134.2 

120.5 

122.4 

6079.567 

c 2 F>z/2 ~ a 2 F>5/2 



2.86 



16.1 

6126.311 

b 2 F 5 / 2 — a 2 F 6 j 2 



11.8 



65.5 

6228.990 

b 2 F5/2 — a 2 Ft /2 



2.03 



10.9 

6148.896 

b 2 F 7 /2 — a 2 F$/2 



3.20 



17.6 

6252.341 

b 2 Fj/2 — Q> 2 Fj/2 



18.8 



100.3 

6652.610 

b 2 Gc)/2 — CL 4j ^7/2 

3.36 

1.97 

2.32 

15.9 

9.3 

11.0 

6727.670 

b 2 G 9/2 - a *F m 

5.12 

4.78 

5.38 

23.7 

22.1 

24.9 

, 9652.716 

b 2 Gq/2 — a 2 Fj/2 

0.71 


1.39 

1.7 


3.3 

9408.356 

b 2 G 9 /2 - a 2 F 5 / 2 

2.62 


1.58 

6.3 


3.8 

6648.980 

b 2 G 7/ 2 - a *Ft /2 

1.87 

1.67 

1.42 

8.9 

7.9 

6.7 

9645.089 

b 2 G 7 /2 — a 2 F 7 /2 

1.02 


0.60 

2.4 


1.4 

7918.452 

b 2 F>z/ 2 ~ a i F 3 /2 



2.62 



8.5 

7977.897 

b 2 F>3/2 - a 4 F 5 / 2 



6.03b 



19.3 

8531.906 

b 2 D 3/ 2 - b 4 F 3/2 



11.00 



30.5 

8587.494 

b 2 D 3/ 2-b 4 F 5/ 2 



5.78 



15.8 

8664.452 

b 2 D'i/2 - b i F 7 / 2 



53.8b 




7979.330 

b 2 D$/2 — CL 4 i ? 7/2 



6.03b 




8568.400 

b 2 D$/2 — b 4 F 7/ 2 



6.44 



17.7 

8663.807 

b 2 D 5 /2 — b 4 F 9/2 



53.80b 



150.2 

8335.590 

6 2 P 3/2 -a 2 F 5/2 



3.74 



10.9 

8526.997 

6 2 P 3/2 -a 2 F 7/2 



13.90 



38.6 

8412.354 

6 2 P 1/2 - a 2 F 5/ 2 



7.72 



22.0 

10069.74b 

b 4 P 5/ 2- ai F 5/ 2 

1.70 


3.57 

3.2 


8.3 


b 4 P 3/ 2-a 4 F 3l 2 







10128.770 

6 4 P 1/2 -a 4 P 3/2 

5.82 


7.84 

13.6 


18.3 


with distance from the central source. Consistently with this 
picture, the novOl slit that was centered on the Sr-filament 
and perpendicular to the radial direction from the central 
source, shows a maximum of the dilution factor at the mid- 
dle of the filament and symmetric drop offs toward both 
sides of this position. 


4 THE TITANIUM AND NICKEL 

ABUNDANCES IN THE SR-FILAMENT 

Once the physical conditions of the Sr-filament have been 
determined we proceed to compute the relative abundances 
of Ni II and Ti II. We use the novOl spectrum because: (1) 
much more spectral interval was covered during that visit 
and the greatly increased number of lines, which permits a 


better statistical determination, and (2) the long aperture 
likely sampled similar regions for the Ni II and Ti II due 
to the tangential slit orientation. Nonetheless, we have con- 
firmed that the computed relative abundances from all other 
observations also agree, within the uncertainties. 

Fig. 10 depicts the N(Ti u)/N(Ni II) as derived from 
the [Ti n] A6648 and [Ni il] A7379 lines. The calculated 
abundance ratio is plotted against the assumed value for 
the dilution factor and for temperatures of 6000 ± 1000 K 
as diagnosed from the spectrum. Given the fact that the 
lines observed from both species have similar excitation en- 
ergies and arise from levels with the same multiplicity and 
parity, they respond similarly to changes in the dilution fac- 
tor. So we see that the uncertainty in the abundance ratio 
is negligible for the mean temperature of 6000 K, and un- 
certainty in the derived abundance does not exceed 20% 
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Figure 5. Ti II line ratios between dipole allowed and forbidden 
line against the logarithm of the dilution factor. The line ratios 
computed at N e = 10®, N e = 10 7 , N e = 10 s cm -3 are presented 
in separate panels located from left to right in this figure. For 
each value of N c the line ratios are computed for temperatures of 
4000, 5000, 6000, and 7000 K. The observed line ratios for each 
observation are indicted by horizontal lines as dashed-line: marOO; 
solid line: novOl 

within the whole range of dilution parameters and temper- 
atures. Thus, we conclude that N(Ti li)/N(Ni il) « 2, with 
a conservative uncertainty that may not exceed a factor of 
three, including the dispersion from all the lines which con- 
tains non-systematic uncertainties in the atomic data and 
the treatment of the pumping continuum field. 

To estimate the total N(Ti)/N(Ni) abundance ratio we 
notice that only neutral and singly ionized ions are present in 
the spectra. It is expected that the ionic ratios Ti i/Ti II and 
Ni i/Ni n are similar since Ti I and Ni I have approximately 
the same ionization potential. Moreover, from observations 
one can conclude that Ni and Ti are both singly ionized, as 
there are no lines from the neutral stages of these ions in the 
spectra. Thus, N(Ti)/N(Ni ) ~ N(Ti + )/N(Ni + ) w 2 in gas 
phase. This value is significantly higher than the best solar 
abundance ratio of ~0.05. In other words, the gas phase 
N(Ti)/N(Ni) abundance ratio in the Sr-filament is nearly 
two orders of magnitude greater than solar value. 

Further, we have checked on the N(Ni n)/N(Fe II) 
abundance ratio in the filament. This calculation is com- 
plicated by the fact that the [Fe IlJ emission comes from a 
more extended region in velocity space. Thus, we only con- 
sider the integrated flux of [Fe il] lines within the range in 
velocity space that coincides the [Ni il] emission we find a 


Figure 6. Ti II line ratios between forbidden lines, /([AA4927 + 
6079] ) //( [A6648] ) , against the logarithm of the dilution factor. 
The line ratios computed at T e = 4000, 5000, 6000, and 7000 K 
and N e = 10 8 , N e = 10', N e — 10 8 cm -3 . The observed line 
ratios for each observation are indicted by horizontal lines as 
dashed-line: marOO; dotted-line: aprOl; solid line: novOl 

value for the N(Ni n)/N(Fe II) ratio roughly consistent 
with the solar N(Ni)/N(Fe) ratio, assuming that the ionic 
fractions of Ni n and Fe II are tied to each other by charge 
exchange and are thus within a factor of two of each other. 
It must be noted, however, that the filament exhibits promi- 
nent emission from Fe I, suggesting a large abundance of this 
ion, but no emission from Ni I or Ti I within the sensitivity 
of the observations. 


5 TITANIUM, NICKEL, AND IRON 
CONDENSATION IN r] Carinae 

The observed large overabundance of titanium with respect 
to both nickel and iron in the gas phase is quite striking. 
Ti/Fe abundance ratios higher than solar by up to 0.4 dex 
are known to occur in metal-poor and extremely metal-poor 
stars (Ryan et al. 1996), but to the best of our knowledge 
there are no reports of such a high overabundance as we 
have analyzed here. It is likely that the abnormal relative 
abundances found here are only representative of the gas 
phase, rather than the total absolute abundances. Under the 
conditions of the Sr-filament nearly 90% of Fe and Ni are 
condensed into grains. Thus, a large Ti/Ni gas phase abun- 
dance ratio would be attained if some kind of destruction 
mechanism of Ti-bearing condensates acted in the filament. 
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Figure 7. Ti II line ratios between forbidden lines vs. T c . 
The line ratios are: /([AA6148 4- 6252])//([A6648]) (Ratio 1), 
I ([AA6652+6727])// ([A6648]) (Ratio 2), /([AA7918+7977+8532+ 
8587])/J([A6648]) (Ratio 3), and /([A100696])//([A6648]) (Ratio 
4). The line ratios computed for lV e = 10 8 , N e = 10 7 , N e = 10 8 
cm -3 are displayed in separated panels from left to right. In each 
case logio iu=-9.0 is indicated by the solid line and log io w=- 8.5 
and -9.5 are indicated by dotted lines. The observed line ratios 
for the novOl are marked with horizontal lines 

Although we do not have a clear explanation for such a 
mechanism we devote this section of the paper to study the 
chemistry of the Ti, Ni, and Fe under the conditions of the 
Sr-filament. 

Condensation temperatures of Ti, Fe, and Ni are some- 
what different and Ti is generally more refractory than Fe 
and Ni during condensation from a solar composition gas 
(e.g., Lodders 2003). As described below, Ti also forms a 
more refractory condensate than Fe and Ni from the overall 
composition relevant for the ejecta here. Condensation tem- 
peratures depend on bulk composition and total pressure. In 
a solar composition gas with C/O = 0.5, titanium condenses 
as calcium titanate (e.g., perovskite, CaTiOa), followed by 
an iron and nickel metal alloy at somewhat lower temper- 
atures (Lodders 2003). Trace elements such as Sr condense 
into solid solution with major element host phases, and con- 
densation of Sr as ’’SrTiOs” into CaTiOa proceeds at tem- 
peratures between that of CaTi 03 and metal alloy forma- 
tion. As long as the overall composition of r] Carinae ’s ejecta 
remained near-solar, solar-type condensates should also have 
formed. However, in the Sr-filament and in the rest of the 
emission ejecta within the Homunculus, the Little Homuncu- 
lus and the Weigelt Blobs, very little oxygen has ever been 


Figure 8. Ti II line ratios between forbidden lines vs. lV e . The line 
ratios are computed for T e = 6000 K (solid line) , and T e = 5000 K 
and T e =7000 K (dashed lines) . The dilution factor was fixed in 
logio W—-9.0. The observed line ratios are indicted by horizontal 
lines as dashed-line: marOO; dotted-line : aprOl; solid line: novOl 


detected. Moreover, in the Weigelt Blobs where [O i] A6300 
emission is seen, oxygen seems to be underabundant by at 
least two orders of magnitude (e.g. Thackeray 1967, Zeth- 
son 2001, Verner et al. 2005). In the outer ejecta carbon and 
oxygen are severely depleted relative to solar abundances 
(Davidson et al. 1982; Dufour et al 1997; Smith and Morse 
2004). Hydrogen is also depleted whereas N and He are en- 
riched. In order to explore the condensation chemistry, we 
adopt a solar composition gas in which solar C and O abun- 
dances are both reduced by a factor of 100. The N abundance 
is increased by the amount removed from the solar C and 
O abundances so that the sum of solar CNO nuclei is pre- 
served. The He abundance is taken as solar plus a quarter of 
the solar H abundance, and the H abundance is set to zero. 
Note that the chemistry of Ti, Fe, Ni, and Si is relatively 
insensitive to the total H abundance because there are no 
major H-bearing gases for those elements. Even if H were 
fairly abundant, the following condensation results would 
change little. Descriptions of the thermochemical compu- 
tation methods are given in Lodders (2003) and Lodders 
& Fegley (1997). The results of this calculation serve as a 
good illustration of the dust chemistry, and should remain 
qualitatively correct despite variations of the precise relative 
abundances. The condensates (Fig. 11) for this ”CNO” pro- 
cessed composition are first TiC, and at about 100 K lower 
temperature, an iron-nickel silicide alloy, (Fe,Ni)Si. Below 
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MarOQ NovOI 



Angular distance (arsec) 

Figure 9. [Ti tl] flux line ratios, /(AA6652 + 6727)//(A4918) (up- 
per plots) and derived dilution factors (lower plots) against the 
angular distance along the spectral slits of marOO (radial across 
the filament from the central source) and novOX (tangential across 
the filament). 


total pressures of ~ 10 -7 bar, SiC condensation tempera- 
tures are within a few degrees of those of the silicide alloy. 

The condensates from the ”CNO” processed composi- 
tion are the same as appear when C and O abundances of an 
otherwise solar composition gas are only slightly (i.e. orders 
of factor 2) modified so that C/O > 1 (Lodders & Fegley 
1995, 1999). Here the C/O ratio remains solar but the large 
decrease in absolute C and O abundances is responsible for 
the appearance of the carbides and silicide because there is 
no oxygen for SiO and silicate or oxide formation. 

The elemental composition here has enough C for con- 
densation of all Ti as TiC (the total Ti/C is ~ 0.03) but 
only a little SiC can form because Si is more abundant than 
C (Si/C ~ 14) and Si condenses also into (Fe,Ni)Si. The 
overall C abundance is so low that graphite does not con- 
dense at the relatively high, temperatures where TiC and 
(Fe,Ni)Si form. If carbon were completely absent here, ti- 
tanium nitride (TiN, osbornite) would condense instead of 
TiC at about the same temperature, SiC would be com- 
pletely absent, and condensation of (Fe,Ni)Si would remain 
unchanged. 

Condensation temperatures describe when a compound 
starts to form (i.e., where the saturation ratio reaches unity) 
but lower temperatures are needed to fully remove an ele- 
ment from the gas into a condensate, as shown for a constant 
total pressure in Fig. 12. Complete titanium removal from 



Figure 10. N {Ti + ) / N{Ni + ) abundance ratio from the spectrum 
of the Sr-fiiament as obtained for any given dilution factor be- 
tween 10~ u and 10 -9 . The solid line depicts the abundance ra- 
tios for an adopted temperature of 6000 K, while the results for 
temperatures of 5000 and 7000 K are indicated by the dotted 
lines. The vertical dashed lines indicate the range of the dilution 
factor as diagnosed from the spectrum 



Figure 11. Condensation temperatures as a function of total 
pressure 
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Table 7. Comparison between observed and theoretical line in- 
tensities for the novOl spectrum. Line intensities are relative to 
/([A6648]. The theoretical spectrum was computed for login = 
-9.0, T e = 6000 K, and N e = 10 7 cm" 3 


A 

/ (A)// ( [A6648] ) 
Obs. Theo. 

4918.26 

15.7 

6.38 

4927.29 

18.2 

18.3 

4984.20 

2.6 

1.01 

6079.57 

2.4 

0.91 

6126.31 

9.7 

7.23 

6148.90 

2.6 

1.78 

6153.62 

3.1 

0.25 

6228.99 

1.6 

0.91 

6252.34 

14.9 

11.70 

6648.98 

1.0 

1.0 

6652.61 

1.6 

1.49 

6727.67 

3.7 

3.72 

7918.45 

1.3 

0.51 

7977.90 

2.9 

0.84 

8335.59 

1.6 

0.80 

8412.35 

3.3 

2.46 

8527.00 

5.7 

4.50 

8531.91 

4.5 

4.39 

8568.40 

2.6 

8.98 

8587.49 

2.4 

8.04 

9645.09 

0.2 

0.98 

9652.72 

0.5 

0.74 

10069 b 

2.2 

2.12 

10128.77 

2.7 

3.26 


the gas requires a drop of about 100 K below the condensa- 
tion temperature of TiC and all Ti is condensed at temper- 
atures when (Fe,Ni)Si condensation begins to proceed. 

Given the nature of the Ti-, Ni-, and Fe-bearings de- 
scribed here we suggest that selective photoevaporation of 
Ti-bearings is a possible explanation for the gas phase Ti/Ni 
abundance found in the filament. Two different scenarios 
may be suggested: (1) there was some spatial separation be- 
tween the TiC and(Ni,Fe)Si condensates, related to their 
different condensation temperatures, that left the TiC rich 
gas more exposed to photoevaporation; (2) the TiC grains 
are generally smaller than the (Ni,Fe)Si grains, and so the 
TiC condensates were much less resistant to photoevapora- 
tion. 

These scenarios sketched out here seem promising to 
explain the observations and need to be studied in more 
detail. Further, it is important to determine the abundance 
ratios of other elements of different volatility. 


6 SUMMARY AND CONCLUSIONS 

We have studied the Ni II and Ti II spectra of the Sr-filament 
in rj Carinae , as taken in three different epochs with the 
HST. In doing this we build spectral models for the two 
ions using the best published data for Ni II, while we carry 
out extensive calculations of radiative and electron impact 
rates for Ti II. 

The spectral models were employed for diagnosing 
the physical conditions of the emitting region. The results 
are generally consistent with previous determinations using 



temperature, K 


Figure 12. The fraction of an element condensed as a function 
of temperature at a constant total pressure of 10~ 7 bar 

[Sr ii] and Sr ii lines. Moreover, the [Ti ii] lines allow us to 
constrain the conditions for observations with the aperture 
across the r/ Carinae ’s central source to N e ~ 10' cm -3 , 
T e — 6000± 1000 K, and dilution factor logiow = -9.0±0.5 
for an assumed blackbody continuum radiation with tem- 
perature of 35 000 K. For observations with the slit perpen- 
dicular to the radial direction the dilution factor seems to 
be. somewhat lower. 

Study of the observed [Ti ii] emission reveals a very 
large overabundance of Ti in gas phase with respect to Ni 
and Fe relative to the cosmic values. 

We study the chemistry of Ti, Ni, and Fe in CNO-cycled 
gas and find that Ti condenses as TiC, or alternatively as 
TiN, while Ni and Fe are locked into (Ni,Fe)Si. Moreover, 
the TiC condenses at higher temperatures than Ni- and Fe- 
bearings. Later, we suggest that the TiC would be selectively 
destroyed by the stellar UV radiation leading to a relative 
overabundance of Ti in gas phase. 

This explanation for the unusual Ti/Ni ratios in the 
Sr-filament requires that significant dust and molecules had 
formed in this material shortly after the Great Eruption, but 
have since been largely photoevporated here. In this context, 
it is very interesting that near-IR emission from molecular 
hydrogen is seen throughout the polar lobes of the Homuncu- 
lus Nebula (Smith 2002a), but this same H 2 emission is not 
seen at comparable levels from the equatorial skirt in the 
vicinity of the Sr-filament (Smith 2002a, 2004). Today, the 
gas densities we infer in the Sr filament (~10 7 cm -3 ) are very 
high compared to many other circumstellar nebulae around 
massive stars, and are similar to the extremely high densities 
seen in the walls of the Homunculus (e.g., Gull et al. 2005; 
Nielsen et al. 2005). A weaker radiation field at larger radii 
is not the only explanation for this difference: H 2 emission 
is seen both in the polar caps of the Homunculus at larger 
radii from the star than the Sr-filament, and in the side walls 
of the polar lobes near the equator, even closer to the star 
than the Sr filament (Smith 2004). In the past when the Sr- 
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filament ejecta were even denser, it is. conceivable that H 2 or 
other molecules may have existed, but have since been de- 
stroyed. If true, this would seem to support the scenario of 
unusual chemistry resulting from selective photoevaporation 
and CNO processed ejecta, as we have suggested here. 
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